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I.      INTRODUCTION 

Finite  state  space  semi-Markov  models  find  applications  in  a  variety  of  fields  such 
as  queueing  theory,  reliability,  and  clinical  trials.  Often  interest  in  the  application  of 
these  models  centers  on  the  distribution  of  the  first  passage  time  to  a  state  or  a  set  of 
states  representing,  for  example,  the  lifetime  of  a  system  or  the  end  of  a  busy  period  of 
a  server.  Suppose  that  the  observations  of  the  path  of  the  semi-Markov  process  are  all 
that  is  known  about  the  process.  The  problem  is  to  estimate  the  probability  that  the  first 
passage  time  has  not  occurred  before  time  t. 

Censored  data  problems  arise  frequently  in  medical,  and  also  in  engineering  system 
reliability  applications.  For  example,  in  medical  survivorship  studies  some  subjects  may 
be  lost  to  follow-up,  or  available  data  may  be  analyzed  before  all  subjects  have  expired. 
In  the  equipment  reliability  context,  observed  units  may  still  be  in  operation,  perhaps 
after  several  previous  failures,  at  the  time  of  the  analysis. 

Three  possible  estimators  will  be  considered.  The  three  estimators  use  different 
amounts  of  information  concerning  the  process.  One  estimator  uses  only  the  observed 
first  passage  times.  Another  estimator  makes  parametric  assumptions  concerning  the 
sojourn  time  distribution  and  uses  maximum  likelihood.  A  third  approach  uses  an  ex- 
ponential approximation  to  the  probability  and  empirical  distributions  to  estimate  the 
sojourn  time  distributions. 

The  three  estimators  were  investigated  for  a  specific  semi- Markov  process  with 
uncensored  data  in  KimfRef.  1]  and  Jacobs[Ref.  2].  Results  of  a  simulation  study  of  the 
three  estimators  using  censored  data  are  reported  in  Gallagher  [Ref.  3].  The  emphasis  in 
this  latter  study  is  on  the  behavior  of  the  point  estimates  with  mean  biases  and  standard 
errors  being  given. 

In  this  thesis  the  investigation  of  the  behavior  of  the  three  estimators  with  censored 
data  is  continued.  The  emphasis  here  is  on  confidence  intervals  for  the  three  estimators. 
In  Chapter  2,  the  three  estimators  are  described  and  the  respective  confidence  interval 
procedures  considered  introduced.  Chapter  3  contains  the  details  of  the  simulation  ex- 
periment and  its  results.  Finally,  conclusions  from  the  study  are  given  in  Chapter  4. 


II.     NATURE  OF  PROBLEM 

A.  PROBLEM 

The  semi-Markov  process  model  considered  is  as  follows.  Suppose  we  observe  N 
individuals.  Let  X,{i)  be  the  state  of  the  i*  individual  at  time  /.  We  will  assume 
{X,(i);t>0}  i  =  1,2,3,. ...,iV.  are  independent  identically  distributed  semi-Markov  proc- 
esses having  the  same  probability  law  as  {X,;t>0}.  The  process  {A^r^O}  is  a  semi- 
Markov  process  with  three  states  {0,  1,2}.  The  individuals  start  at  t  =  0  in  state  1. 
Upon  leaving  state  1,  the  process  transitions  to  state  0  with  probability  6  and  to  state  2 
with  probability  1  —  9.  From  state  2  the  process  transitions  to  state  1  with  probability 
1.  State  0  is  an  absorbing  state.  The  first  passage  time  to  state  0  will  be  referred  to  as 
the  time  of  death.  The  N  individuals  are  censored  independently.  The  censoring  times 
are  exponentially  distributed  with  a  mean  of  — .  The  entire  path  of  transitions  and 
sojourn  times  are  observed  until  the  time  of  censoring  or  death. 

Let 

D  =  inf{t>0;Xt  =  0} 

and 

S(t)  =  P{D>t) 

where  D  is  the  time  of  death  (or  entrance  to  state  0).    The  problem  is  to  estimate  the 
survival  probability  P{D  >  t)  with  the  censored  data  of  N  individuals. 

B.  ESTIMATORS 

1.     Kaplan-Meier  Estimator 

A  non-parametric  estimator  of  the  distribution  function  for  censored  data  is  the 
Kaplan-Meier  estimator(K.M.E.)  which  is  often  called  the  product  limit  estimator  [Ref. 
4].  Let  L\,  U2,  —,Un  be  independent  identically  distributed  random  variables  with  distrib- 
ution G  having  a  density  function.  Let  Vu  V2, ...,  Vn  be  independent  identically  distrib- 
uted times  to  censoring  with  a  continuous  distribution  function.  Let 

Z,  =    min(I/„  V,), 
and 


JO  ifU^V, 

[I  Otherwise. 

The  Z,  are  the  observed  times  and  <5,  is  an  indicator  of  whether  or  not  the  i'h  observation 
is  censored.  Let  Z(1)  <  ZQ)  <  ...  <  Zw  be  the  order  statistics  of  {Z,}  and  <5(0  be  the  corre- 
sponding values  of  {<5,}.  It  is  assumed  that  there  will  be  no  ties  since  the  underlying 
distribution  functions  are  continuous.  The  Kaplan-Meier  estimate  of  the  survival  func- 
tion S(t)  =  (l-G(r))  is 


n 


(n  -  i) 


l-<5<0 


if    t<Z{n) 


(«-/+  1) 

S(t)     =     <!     [/:Zw<r]  (2.1) 

I  0  if    t>Z{n)&6{n)  =  0 

Undefined  if    t>  Z^  &  S^  =  1 . 

A 

The  variance  of  S(t)  is  given  approximately  by 


A       ,2 


y s.l- 

/_j    {n-  i){n  -  i 


Varmn^lSW     )      („_l)(„_l  +  l)  (2-2) 

[Ref.  4:  p.  464].  The  Kaplan-Meier  estimator  using  the  death  times  of  the  N  individuals 
will  be  denoted  by  PK{D  >  /}. 

If  the  Kaplan-Meier  estimate  is  undefined,  we  investigate  the  effect  of  two 
methods,  defined  as  MOD  1  and  MOD  2,  to  make  the  Kaplan-Meier  estimate  honest. 
MOD  1  and  MOD  2  are  defined  as  follow.  In  MOD  1,  the  remaining  mass  of  the  esti- 
mated survival  function  is  assigned  to  the  last  datum  Z(n)( which  is  censored).  In  MOD 
2,  if  the  last  k  data  points  are  censored  the  remaining  mass  of  the  survival  function  is 
distributed  equally  among  the  k  data  points.  For  example,  if  the  estimated  survival 

probability  at  the  last  uncensored  point  is  0.2  and  there  are  additional  two  data  points 

0  2 
which  are  censored,  then  each  of  these  additional  points  is  assigned  a  mass  of  — r—  . 


2.     Maximum  Likelihood  Estimator 

In  this  subsection,  the  maximum  likelihood  estimator(M.L.E.)  will  be  given  for 
the  special  case  in  which  the  sojourn  time  in  state  /  is  exponentially  distributed  with 
mean  -£-(/=  1,  2). 

Let  R,j  be  the  number  of  transitions  from  state  i  to  state  j  for  one  individual. 
The  log  likelihood  function  for  an  individual  is 

/  =  Rn  ln(l  -  0)  +  Rw  In  6  +  R2l  In  p2  +  {Rl0  +  Rn)  In  Pl  -pjx-  p2T2        (2.3) 

where  T,  (/  =  1,2)  is  the  total  time  spent  in  state  i  before  death  or  censoring. 

The  maximum  likelihood  estimators  using  the  data  from  all  N  individuals  are 


where 


and 


RXQ  +  Rn 


g-  ~  "10~    ;  (2.4) 


A  /Cin  +  An 

Pi"      ~       ;  •  (2.5) 


i 


p2  =  ^-;  (2-6) 


A       #  death  times  that  are  censored 

c=- - - ,  (2.7) 

T,  +  T, 


N 

Rq-J^Rfo)  (2.8) 


N 

T^Y.Un)  (2.9) 

with  R0(n)  being  the  number  of  transitions  from  i  to  j  for  the  n'H  individual  and  T{n)  being 
the  total  time  spent  in  state  i  before  death  or  censoring  for  the  i*  individual. 
To  obtain  asymptotic  variances  for  these  estimators,  note  that 


-c2l 

#10 

K 

12 

ee2 

02 

+ 

(1- 

o? 

-c2l 

#10 

+  #12 
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Pi 

—  0  I 

7?21 

a    2 
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2 
P2 

Further 


n=\ 


=  Y (i  -  d)    p!_    [(    p!_    )( — ^—  )(i  -  e)f 

L-T  '  Pj  +  c   LV  pt  +  c  /v  p2  +  c 

M«l 

(l-0)p,(p2  +  c) 


0PiP2  +  c[pj  +  p2]  +  c2 


rt=1 


=y[(1_e)_4 — a_j. 

/     /  Pi   +  C      P2  +  C    J 

(l~g)PlP2 

0PlP2  +  cl>l  +P2]  +  c2 


^PlP2  +  cCPl  +  P2]  +  C2 


(2.10) 


(2-11) 


(2.12) 


£[*i2]  =  2/</^«}  (2-13) 


£I#2l]=X/5(/?21^W  (2-14) 


£[/?  ]  =  V    p'     0[(i  -  0){  - p-^r r }]"  (2-15) 

p,0(p2  +  c) 


where  c_1  is  the  mean  of  the  exponential  censoring  time.    Thus  for  N  individuals  the 
asymptotic  variances  of  the  estimators  are 

Var${\  =  I{px)-\  (2.16) 

Varfal  =  I{p2YX 
where 

D  D 

];  (2.17) 


(2.18) 
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■  N 
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El 
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Pi 
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El 

*21 
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Pi 
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(2.19) 


The   expression  for  the   survival  function   S{t)  =  P{D  >  t)   for  this   continuous  time 
Markov  chain  is 

s(0  .  {  h±lL  eh<  _  hllJL  eh<)e  _^L_  (2.20) 

/. 2  >■!  /,  -  >.2 

where  /,  X2  are  the  roots  of  the  equation 

dplp2+y(pl+p2)+y2  =  0.  (2.21) 

A 

The  maximum  likelihood  estimator,  denoted  as  PM{D>t),  for  the  survival 
probability  is  [Ref.  5:  p.  5  eqn  1.17] 

PM{D>^Jhr\H^^-H^H  (2.22) 

/^  -  >-2     t         >-2  ^1  J 

where  Xx  and  A2  are  roots  of  the  equation 

eKpi  +  ><pi  +  p2)  +  y2  =  o.  (2.23) 


Since  the  maximum  likelihood  estimators  are  orthogonal,  the  asymptotic  variance  of 
PM{D  >  t]  [Ref.  5:  p.  5  eqn  1.19]  is  approximately 

Var\_PM{D>t)\d,p„Pl\ 

=  Var{d){^T?  +   VarCp,)(-^-)2  +   Var{p2){-^-f  (2.24) 

dd  dp,  dp2 

where 
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+  TT~  "ir-  +  "7T~  "7T~  '  (2-26) 


cpj        Pi        o/j    8px        o).2    dpx 


dS            Op]       r    1      v       1      Vij    dS    d>-\        dS    o?.2 
— —  =  —, — -r-r-  e  2  —  "T—c  '  J  +— : h  — : — ;  (2.27) 

dp2     (;.,->.2)    >.2         >.,  a;.!  ap2     a>.2  aP2 


as  -S        .        9p\       r  i°2     v      ,  P2    .   n.  ;.,rn  010^ 


-^-  = +       °PX       [  -  —  e*  +  ( ■£-  +  Dfe^l-  (2  29) 

a;.2     (;.,  -  ;.2)  +  (;.,  -  ;.2)  L     A2      + i  ;.2  +  ^  J.  ^-^ 


^•/         -  P1P2 


a#     Pi  +  p2  +  2;.,- ' 
a;.,       -xl-ePl 


dPi       P\+  P2  +  2//  ' 

a;.,  _   -^-g/p, 

<?P2       Pi  +  P2  +  2;-/  ' 


(2.30) 


(2.31) 


(2.32) 


3.     Asymptotic  Renewal  Estimator 

In  this  subsection,  we  describe  an  asymptotic  renewal  estimator(A.R.E.)    for 
P{D  >  t}.   A  conditioning  argument  yields  the  following  equation  for  P{D  >  t]  ; 


P{D  >t)=  g(i)  +  (1  -  0)J  (F,*Fa)  (<fc)  P{Z)  >  t  -  s]  (2.33) 


with 


g(r)  =  fi(/)  +  (1  -  0)  [Vj  (A)  F2(r  -  5).  (2.34) 

o 

Thus,  P{Z)  >  t]  satisfies  a  renewal-type  equation  with  defective  inter-renewal  distribution 
L{t)  =  (l-9)(F,*F2)(t)  (2.35) 

where  Ft  is  the  sojourn  time  distribution  in  state  i,  F,(t)  =  1  —  F,(t),  and  {F*F2)(t)  denotes 
the  convolution  of/7!  and  F2.    Following  Feller  [Ref  6],  let  k  be  such  that 


ekl  L(dt)  =  (1  -  0)  <j>x{k)  (j)2{k)  =  1  (2.36) 


where 


0KO-J     eltFidt).  (2.37) 


o 
Then,  under  certain  integrability  conditions,  if  (F1*F2)  is  not  arithmetic 

lim  ekt  P{D>t}=-j-  (2.38) 

where 


ix  =  I     ie*1  I(<fc)  (2.39) 


and 


ks 


b=        eK'g(s)ds 


(2.40) 


ks  r? 


[/fa) +  (1-0)  J     Fx{du)F2{s-u)1ds 


J  [{*,(*)  -  1}  +  (i  -  0)  W*)W)  "  1}]. 


Since 


(i-W,(*)*2(*)-i 


it  follows  that 


b  =  T^-V  +  JJi^^-V 


=  |[^)-i]  +  |[i-^] 


=  y[0iW-l]  +  y  [1-(1 -*)*,(*)] 


=  y  [*,(*)  -  1]  +  y  [1  -  *,(*)  +  e  0,(*)] 


=f<t>m 


(2.41) 


Let  F,  be  the  Kaplan-Meier  estimate  of  F„  and  0  be  the  maximum  likelihood 
estimate  of  6  [Ref.  5:  p.  10];  put 


Too 


US)  = 


elsF{ds). 


(2.42) 


The  asymptotic  renewal  estimator  (PA{D>t))  [Ref.  5:  p.  11  eqn  3.11]  of  the  survival 
probability  P{D>  t)  is 

A 

PA{D>t)    -    **'-*-  (2-43) 

where  A:  is  the  solution  to  the  equation 

U-hfaWiW  -  i;  (2-44) 


M    =    (1-0) 


A  A 


ekss(F^F2)(ds);  (2.45) 


^o 


A  AAA 

b    =    -?-0i(*).  (2-46) 


In  the  simulation  A:  is  obtained  by  using  the  golden  section  search  method. 

The  asymptotic  renewal  estimator  is  undefined  if  all  the  sojourn  times  for  a 
particular  state  are  censored  since  the  Kaplan-Meier  estimator  is  not  defined  in  this  case. 

C.     CONFIDENCE  INTERVAL  PROCEDURES 

A  confidence  interval  for  an  unknown  parameter  gives  both  an  indication  of  the 
numerical  value  of  the  unknown  parameter  and  a  measure  of  how  confident  we  are  of 
that  numerical  value.  Two  statistics  L  and  U  form  a  (1  —  a)  100%  confidence  interval  for 
6,  if  under  repeated  random  sampling  L  <  0  <  U  (1  —  a)  100%  of  the  time.  A  confi- 
dence interval  procedure  for  an  unknown  parameter  0  is  also  used  to  make  a  decision 
concerning  0,  as  in  classical  hypothesis  testing  or  decision  making,  or  to  indicate  the 
accuracy  and  variability  of  a  point  estimator  0. 

Confidence  interval  procedures  for  the  three  estimators  for  P{D  >  /}  will  be  de- 
scribed in  this  section,  starting  with  the  confidence  interval  procedure  for  the  Kaplan- 
Meier  estimator(K.M.E.).  Procedures  for  the  maximum  likelihood  estimator(M.L.E.) 
and  the  asymptotic  renewal  estimator  will  then  be  discussed. 

A  preliminary  transformation  to  approximately  symmetrize  the  sampling  distrib- 
ution of  the  estimator  is  often  beneficial  [Ref.  7].  For  this  study  we  consider  two  tran- 
sformations, the  arc-sine  and  log  transformations,  which  tend  to  stabilize  and  also 
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approximately  symmetrize  the  data.  These  transformations  were  suggested  by  the  work 
of  Gaver  and  Miller  [Ref.  8].  Since  the  individuals  are  independent,  the  number  of  in- 
dividuals surviving  a  fixed  time  t  would  have  a  binomial  distribution  if  there  were  no 
censoring.  The  logarithmic  and  arc-sine  transformations  have  been  beneficial  in  this 
case. 

The  confidence  intervals  for  the  Kaplan-Meier  estimator  and  the  maximum  likeli- 
hood estimator  are  the  asymptotic  normal  confidence  intervals  using  the  transformed 
estimator. 

The  confidence  interval  for  the  arc-sine  transformed  estimator  is  computed  as  fol- 
lows. Let  S(t)  be  the  estimator  of  S(t)  =  P{D  >  t)  .  Since 

d     •  -l   / —  1    1  1  ,~  *-,-. 

—-sin    V*     =    T     } -7=^,  (2-47) 

dx  2    J\  -  x     Jx 

a  Taylor  expansion  yields 


sin%/S(f) )  *    sin-'cV^J")    +    T     /  =- (S(i)  -  S(t)).  (2.48) 

x  2  V(i-s(/))S(/) 


Hence,  the  approximate  variance  of  sha~l(yjS(t) )  is 


Vari  sm\y/I(t)  )]  *    -J       ,    l     ,       VarlS(t)l  (2.49) 

4    (l-S(f))S(/) 


If  S(t)    =    1,  we  set  Var[  sir\-l(J  S(i)  )]  equal  to  0.  Therefore,  a  (1  -  a)  100%  confidence 
interval  for  sin_1v//5(/)  is 


.n~.  ..     .  _i.  / 


un-\y/S{t))   ±    [zl-aidy/Varlun~\yjS{t))]  (2.50) 

where  zx_%jl  is  the  (1  -  a/2)  100%  point  of  a  standard  normal. 

A  confidence  interval  for  the  log  transformed  estimator  of  S(t)  is  computed  as  fol- 
lows. Since 

-^In*   =    \,  (2.51) 

a  Taylor  expansion  yields 

lnS(t)  *  InS(t)    +    -±^{S(t)-S{t)).  (2.52) 


11 


Thus  the  approximate  variance  of  In  S(t)  is 

VarllnS(t)-]  a    -^—  VafS(*)l  (2-53) 

S(t)2 

If  £(/)  =  1,  we  set  Far[S(/)]  =  0  .   A  (1  -  a)100%  confidence  interval  for  In  S{t)  is 


toS(f)    ±    {z^y/VarllnSit)-]  (2.54) 

where  Zi_,/2  is  the  (1  —  a/2)100%  point  of  a  standard  normal. 

1.  Asymptotic  Normal  Confidence  Intervals  for  K.M.E. 

A 

The  asymptotic  variance  of  the  Kaplan-Meier  estimator  (PK{D>t})  is  given 
by  equation  (2.2).  The  asymptotic  normal  confidence  intervals  for 
sin.'1  J PK{D  >  t}  and  In  PK{D  >  t)  are  evaluated  using  equations  (2.50)  and  (2.54)  re- 
spectively. The  corresponding  confidence  intervals  for  PK{D  >  r}  are  formed  by  inverse 
transformation.  If  the  lower  limit  is  less  than  0  it  is  set  equal  to  0  and  if  the  upper  limit 
is  greater  than  1  it  is  set  equal  to  1. 

2.  Asymptotic  Normal  Confidence  Intervals  for  M.L.E. 

The  asymptotic  variance  for  the  maximum  likelihood  estimator  (PM{D  >  r})  is 
given  by  equation  (2.22).  The  asymptotic  normal  confidence  intervals  for 
sin-l(yj  PM{D  >  t]  )  and  In  PM{D  >  t]  can  be  constructed  using  equations  (2.50)  and 
(2.54)  respectively.  Confidence  intervals  for  PM{D  >  t)  can  be  obtained  by  inverse 
transformation.  If  the  resulting  confidence  interval  has  a  lower  limit  less  than  0  it  is  set 
equal  to  0  and  if  it  has  an  upper  limit  greater  than  1  it  is  set  equal  to  1. 

3.  Jackknife  Procedure  for  A.R.E. 

The  jackknife  technique  was  first  introduced  by  Quenouille  [Ref.  9]  and  later 
utilized  by  Tukey  [Ref.  10]  for  bias  reduction  and  robust  interval  estimation.  A  review 
can  be  found  in  Miller  [Ref.  1 1].  The  jackknife  is  designed  to  do  various  jobs  fairly  well, 
however  it  is  desirable  to  avoid  (in  jackknifing)  sampling  distributions  with  (i)  abrupt 
ends  and  (ii)  one  or  more  straggling  tails,  and  it  is  probably  desirable  to  avoid  those  that 
are  strongly  unsymmetrical  [Ref.  12].  Confidence  intervals  for  the  asymptotic  renewal 
estimator  will  be  obtained  using  the  jackknife  procedure  on  the  arc-sine  and  log  trans- 
formed estimates. 

The  jackknife  procedure  for  the  confidence  interval  of  the  asymptotic  renewal 
estimator  (PA{D  >  t})  is  implemented  as  follows. 
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1.  Generate  data  for  N  individuals. 

A 

2.  Compute  PA{D  >  t}  using  all  data. 

3.  Transform  PA{D  >  t}  into  In  PA{D  >  t}  and  sin-1^ PA{D  >  t}  which  will  be  denoted 

as  FLaI„    and  YStU  respectively.    Divide  the  N  individuals  into  n  subgroups  such 

N 
that  each  subgroup  contains  —  individuals. 

A 

4.  Compute  PA{D>t)  leaving  out  all  data  of  the  /'*  subgroup  and  transform  it  into 


In  PA{D  >  t)  and  sin-1  ^j PA{D  >  t}  which  will  be  denoted  as  yl(,  and  ys,  respectively. 

5.  Compute  the  pseudo-values(denoted  as  FL.„  and  YS., ); 

YL-i   =    nYLall    -    (n-l)ylh  (2.55) 

YS*t    =    nYS,n    -    (n-l)ys,  (2.56) 

6.  Compute  the  average  of  pseudo-values  which  is  the  jackknifed  estimate  for  the 
transformed  asymptotic  renewal  estimator(denoted  as  YL.,   and  YS.  ); 


YL*    =    jrlYL*1    +     ...     +    YL*„l  (2.57) 


YS*    =    -jflYS.,    +     ...     +   rS.J,  (2.58) 


and   compute   the   variance   of  the   average   of  the   pseudo-values(denoted   as 
VU,   and  VSi)\ 


n 


yd  -  i^h)Y}YU'  ~  YLt  (259) 
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ft 


v$*    =     ..,..     ^   /  {YS*t         YS.y.  (2.60) 

7.     Compute  the  (approximate)  two-sided  (1  —  a)  100%  confidence  intervals  for  each 
of  the  transformed  estimators  as  follows; 


YL.     ±     t{x^i2){n-\)sjvd  (2.61) 


and 


YS.     ±     tQ^pfn-Dy/VSi  (2.62) 


where  fj_B/2(«  -  1)  is  the  (1  -  a/2)100%  point  of  student's  t  with  n-1  degree  of  free- 
dom. 

Confidence  intervals  for  P{D>t)  are  obtained  by  inverse  transformation.  If  the 
resulting  interval  has  a  lower  limit  less  than  0  it  is  set  equal  to  0  and  if  it  has  an 
upper  limit  greater  than  1  it  is  set  equal  to  1. 
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III.     ANALYSIS  OF  SIMULATION  RESULTS  FOR  CONFIDENCE 

INTERVAL  PROCEDURES 

A.     SIMULATION 

All  simulations  were  carried  out  on  an  IBM  3033AP  computer  at  the  Naval  Post- 
graduate School  using  the  LLRANDOM  II  random  number  generator  package[Ref. 
13].  Plots  of  simulated  estimates  and  confidence  intervals  were  produced  by  an  exper- 
imental APL  package  GRAFSTAT  which  the  Naval  Postgraduate  School  is  using  under 
a  test  agreement  with  IBM  Watson  Research  Center,  Yorktown,  Heights,  NY. 

The  data  for  the  simulation  experiments  are  generated  as  follows; 

1.  An  individual  starts  in  state  1  at  time  0. 

2.  An  exponential  censoring  time  with  mean  1/c  is  generated. 

3.  An  exponential  sojourn  time  in  state  1   with  mean  1/  px  is  generated. 

4.  The  sojourn  time  and  censoring  time  are  compared;  if  the  sojourn  time  is  smaller, 
then  the  sojourn  time  is  recorded  and  given  an  uncensored  index  '0';  if  the  censor- 
ing time  is  smaller,  the  sojourn  time  truncated  at  the  censoring  time  is  recorded  and 
given  a  censored  index  T.  In  the  latter  case,  the  death  time  is  recorded  as  the 
truncated  sojourn  time  and  associated  with  a  censored  index  of  T  and  the  simu- 
lation for  the  first  individual  is  completed. 

5.  If  the  process  is  not  censored  in  state  1,  a  uniform  random  number  is  generated 
and  compared  with  6\  if  less  than  6  the  process  jumps  to  state  0,  and  the  uncen- 
sored death  time  with  index  '0'  is  recorded  and  the  simulation  for  the  first  individual 
is  completed;  if  greater  than  0,  the  process  jumps  to  state  2. 
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6.  An  exponential  sojourn  time  for  state  2  with  mean  1/  p2  is  generated  and  the  total 
time(sojourn  time  in  state  1  plus  sojourn  time  in  state  2)  is  compared  to  the  cen- 
soring time  with  the  same  actions  as  listed  above. 

7.  If  the  process  is  not  censored  by  the  end  of  the  sojourn  time  in  state  2,  the  process 
jumps  to  state  1  and  continues  until  an  uncensored  or  censored  death  occurs.  The 
time  is  recorded  and  the  next  individual  is  started. 

This  procedure  continues  until  N  individuals'  data  have  been  generated.  Using  this 
data,  the  Kaplan-Meier  estimate  (PK(t))  ,  maximum  likelihood  estimate  (P^t)), 
asymptotic  renewal  estimate  (PA{t)),  and  their  respective  90%  and  80%  two-sided  confi- 
dence intervals  are  computed.  This  completes  the  one  super-replication.The  simulation 
is  replicated  for  SR  =  500  super-replications  utilizing  different  seeds  to  generate  the 
data.  For     the     simulated     model     described     above     parameter     values     of 

Pt  =  1,  p2  =  1,  6  =  0.5,  and  c  =  —  ,    —  ,    -tt-  ,  and  -r-rr-  are  used,  and  for  each  super- 
replication  N  =  50  individuals'  data  are  generated. 

B.      ANALYSIS 

In  this  section  results  from  the  simulation  experiments  will  be  reported.  In  Ap- 
pendix A,  the  true  survival  probability,  which  is  obtained  using  equation  (2.20),  is  given 
at  the  various  values  oft.  The  times  considered  are  t  =  1.0,  3.0,  5.0,  7.0,  9.0,  11.0,  13.0, 
and  15.0. 

Tables  2  through  5  of  Appendix  B  report  the  number  of  super-replications  that  still 
have  defined  Kaplan-Meier  estimates  at  various  times  t.  The  column  headed  by  KM0 
reports  the  number  of  super-replications  that  have  defined  Kaplan-Meier  estimates  at 
time  t  for  the  death  time  survival  function.  The  columns  headed  KM1.  KM2,  repectively 
show  the  number  of  super-replications  still  defined  for  the  Kaplan- Meier  estimate  of  the 
survival  function  at  time  t  for  the  sojourn  times  in  state  1,  respectively  state  2.  These 
numbers  indicate  the  effect  the  censoring  has  on  the  Kaplan-Meier  estimate.  As  ex- 
pected, increasing  the  mean  time  to  censoring  increases  the  number  of  super-replications 
censored.  In  all  cases  the  survival  functions  for  the  sojourn  times  in  state  1  and  2  are  less 
heavily  censored  than  that  of  the  death  times. 

Plots  for  comparing  the  methods  investigated  for  the  Kaplan-Meier  estimator  are 
given  in  Figures  1  through  3.  The  difference  between  the  mean  of  the  estimated  survival 
probability  and  the  true  survival  probability  (S{t)  —  P{D  >  /}),  the  relative  differences 
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({S(t)  -  P{D  >  t})jP{D  >  i}),  and  the  relative  root  mean  square  errors 
(RMSE/P{Z)  >  /})  are  respectively  plotted  in  Figures  1  through  3.  The  numerical  results 
are  recorded  in  Tables  6  through  17.  The  means,  and  root  mean  square  errors(RMSE) 
from  SR  =  500  super- replications  are  computed  as  follows; 


SR 


SR 


D<( 


s(t)  -   -±r  )  s,(D 


and 


[~SR~}  (5/(0  "  5(/))2]l" 


where  S,(t)  is  the  point  estimate  of  the  true  value  S(t)  at  time  t  in  the  /"'  super-replication 
and  SR  is  the  number  of  super-replication. 

If  a  Kaplan-Meier  estimate  is  undefined,  its  value  is  taken  to  be  its  last  defined 
value.  The  modifications  to  the  Kaplan-Meier  estimate,  MOD  1  and  MOD  2,  are  de- 
scribed in  Chapter  2.  Since  the  unmodified  Kaplan-Meier  survival  function  may  not 
tend  to  0  as  t  ->  oo  ,  it  tends  to  overestimate  the  true  survival  probability.  On  the  other 
hand,  since  the  modifications  MOD  1  and  MOD  2  make  the  undefined  Kaplan-Meier 
survival  function  honest,  they  do  well  for  large  t  and  have  less  variability  than  the  un- 
modified K.M.E..  However  the  modifications  appear  to  bias  the  estimates  for  moderate 
times  t.  All  these  methods  improve  as  the  mean  censoring  time  decreases.  MOD  1  does 
slightly  better  than  MOD  2. 

In  Figures  4  through  9,  plots  are  presented  for  the  comparison  between  the  three 
estimators.  The  differences,  relative  differences,  and  relative  RMSE's  are  computed  and 
plotted  in  the  same  manner  as  before.  Figures  4  through  6  show  results  for  the 
Kaplan-Meier  and  asymptotic  renewal  estimators  using  the  unmodified  Kaplan-Meier 
estimator.  Figures  7  through  9  show  the  results  obtained  when  the  Kaplan-Meier  esti- 
mates are  modified  to  be  MOD  1  in  both  the  K.M.E.  and  A.R.E..  The  numerical  results 
are  recorded  in  Appendix  C.  As  expected,  the  maximum  likelihood  estimator(M.L.E.), 
which  uses  the  most  correct  information  about  the  model,  tends  to  have  the  smallest 
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relative  RMSE's  and  differences  between  estimated  and  true  means.     In  the  case  of 
greatest  censoring,  the  M.L.E.  shows  a  slight  bias. 

The  Kaplan-Meier  estimator(K.M.E.),  which  uses  the  least  information  about  the 
model,  does  well  for  small  times  even  in  the  greatest  censoring  case,  but  it  is  worse  than 
the  others  for  large  times  due  to  the  undefined  estimates.  It  has  the  largest  values  of 
mean  square  error  among  the  three  estimators.  Therefore  the  K.M.E.  tends  to  be  more 
variable  than  the  other  estimators. 

Not  surprisingly,  the  asymptotic  renewal  estimator(A.R.E.),  which  uses  an 
asymptotic  exponential  approximation  to  the  probability,  has  high  relative  RMSE's  and 
differences  between  estimated  and  true  means  for  small  times  t.  As  t  increases  the  rela- 
tive RMSE's  and  differences  between  the  means  decrease.  Tables  6  through  9  indicate 
that  for  moderate  to  large  times  the  A.R.E.  has  RMSE's  that  are  comparable  or  less 
than  those  of  the  M.L.E..  The  time  at  which  they  become  comparable  is  a  function  of 
the  amount  of  censoring.  If  the  mean  censoring  time  is  larger,  the  RMSE's  become 
comparable  sooner.  Tables  6  through  12  indicate  that  while  modifying  the  Kaplan-Meier 
estimate  of  the  survival  function  of  the  death  times  improves  it  for  large  times  t,  it  cre- 
ates a  bias  for  moderate  times  t.  Using  modified  Kaplan-Meier  estimators  of  the  survival 
functions  for  the  sojourn  times  in  states  1  and  2  does  not  improve  the  asymptotic  re- 
newal estimator. 

In  Figures  10  through  17,  plots  are  presented  for  comparing  the  confidence  interval 
procedures  of  the  three  estimators.  In  order  to  compare  the  performance  of  the  confi- 
dence interval  procedures,  we  use  the  following  measures:  the  coverage  fraction(CVR), 
the  average  half  length(AHL)  of  the  confidence  interval,  and  the  standard 
deviation(SHL)  of  the  half  lengths.  The  AHL  is  determined  by  summing  the  half  lengths 
of  the  confidence  intervals  of  the  replications  and  dividing  that  sum  by  the  number  of 
replications(SR).  The  SHL  is  computed  by  summing  the  square  of  the  differences  be- 
tween the  individual  half  lengths  and  the  average  half  length,  dividing  that  value  by 
(SR-1)  and  taking  the  square  root  of  the  resulting  value.  Among  confidence  interval 
procedures  which  achieve  the  desired  coverage  rate(0.90,  and  0.80),  the  confidence  in- 
terval procedure  which  yields  the  smallest  AHL  is  preferred.  Also  preferred  is  a  small 
SHL  representing  a  stable  confidence  interval  procedure. 

For  each  procedure,  the  number  of  intervals  covering  the  true  value  P{D  >  t}  is 
recorded  as  well  as  the  number  of  intervals  that  are  too  high  or  too  low.  These  results 
are  reported  in  Tables  18,  20,  22,  24,  26,  28,  30,  and  32  (Appendix  D).  Next  to  each 
coverage  count  is  given  the  corresponding  coverage  fraction  in  parentheses.      If  a 
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(1  —  a)  100%  confidence  interval  procedure  is  performing  well,  then  this  interval  should 
cover  about  (1  —  a)100%  of  time.  A  95%  confidence  interval  for  the  coverage  fraction 
is  computed  using  P  ±  1.96[P(1  —  P)ISR]J  ,where  P  is  the  proportion  of  (1  —  a)%  con- 
fidence intervals  that  cover  the  true  value  of  P{D  >  t)  and  SR  is  the  number  of  super- 
replications.  Therefore,  if  a  80%(90%)  confidence  interval  procedure  is  working  well, 
then  out  of  500  super- replications  between  382(436)  and  418(463)  confidence  intervals 
should  cover  the  true  value.  These  intervals  correspond  to  0.7649(0.8351)  and 
0.8737(0.9263)  respectively.  The  average  half  lengths  of  the  confidence  intervals  for 
S(t)  are  reported  in  tables  19,  21,  23,  25,  27,  29,  31,  and  33  (Appendix  D).  The  standard 
deviation  of  the  half  length  is  given  in  parentheses  below  the  average  half  length.  If  an 
estimator  is  performing  well,  its  confidence  interval  should  not  only  have  the  correct 
coverage  fraction  but  also  a  small  average  half  length. 

In  Figures  10  through  13,  the  coverage  fraction  and  the  average  half  length  of 
confidence  intervals  using  the  log  transformation  for  each  estimator  are  presented.  In 
Figures  14  through  17,  the  coverage  fraction  and  the  average  half  length  of  confidence 
intervals  using  the  arc-sine  transformation  are  presented.  The  two  horizontal  lines  in  the 
coverage  plot  show  the  95%  confidence  interval  for  the  coverage  fraction. 

The  following  remarks  concern  the  confidence  intervals  obtained  using  the  arc-sine 
transformation.  The  asymptotic  normal  confidence  intervals  for  the  M.L.E.  have  the 
correct  coverage  for  small  to  moderate  times  but  tend  to  slightly  undercover  for  large 
times.  The  asymptotic  nomal  confidence  intervals  for  the  K.M.E.  have  the  correct  cov- 
erage for  small  times  but  undercover  for  moderate  to  large  times  due  to  the  undefined 
estimates.  As  the  mean  censoring  time  increases,  the  coverage  of  the  Kaplan-Meier 
confidence  intervals  improves.  The  confidence  intervals  for  the  asymptotic  renewal  esti- 
mate tend  to  undercover.  Using  a  modified  Kaplan- Meier  estimate  makes  very  little  dif- 
ference in  the  results. 

The  following  remarks  concern  the  confidence  intervals  obtained  by  using  the  log 
transformation.  Once  again  the  confidence  intervals  for  the  K.M.E.  undercover  for 
moderate  and  large  times  t.  The  coverage  is  slightly  worse  than  that  obtained  using  the 
arc-sine  transformation.  The  confidence  intervals  for  the  M.L.E.  have  the  correct  cov- 
erage and  the  smallest  AHL.  The  confidence  intervals  for  the  Jackknifed  A.R.E.  have 
the  correct  coverage  for  all  but  the  80%  interval  at  t=  1  for  c=  1/2.  The  average  half 
lengths  of  the  A.R.E.  confidence  intervals  are  larger  than  those  of  the  M.L.E.  confidence 
intervals.  Using  a  modified  Kaplan-Meier  estimate  makes  very  little  difference  in  the  re- 
sults. 
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Figure  1.        COMPARISON  BETWEEN  METHODS  FOR  K.M.E.:      DilTerence 
between  Mean  and  True  value. 
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Figure  2.        COMPARISON  BETWEEN  METHODS  FOR  K.M.E.:       Relative 
DnTerence  between  Mean  and  True  value. 
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Figure  3.         COMPARISON  BETWEEN  METHODS  FOR  K.M.E.:       Relative 
Root  Mean  Square  Error. 
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Figure  4.       COMPARISON  BETWEEN  ESTIMATORS(UNMOD):      Difference 
between  Mean  and  True  value. 
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Figure  5.         COMPARISON  BETWEEN  ESTIMATORS(UNMOD):       Relative 
Difference  between  Mean  and  True  value. 
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Figure  6.        COMPARISON  BETWEEN  ESTIMATORS(UNMOD):       Relative 
Root  Mean  Square  Error. 
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Figure  7.        COMPARISON  BETWEEN  ESTIMATORS(MOD  1): 

between  Mean  and  True  value. 
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Figure   10. 


CONFIDENCE   INTERVAL(90%,  LOG,  UNMOD): 

Fraction  and  Average  Half  Lensth. 
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Figure  11. 
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Figure  12.  CONFIDENCE   INTERVAL(90%,   LOG,   MOD    1): 

Fraction  and  Average  Half  Lensth. 
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CONFIDENCE   INTERVAL(80%,   LOG,   MOD    1): 

Fraction  and  Average  Half  Length. 
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Figure  14. 


CONFIDENCE  INTERVAL(90%,  ARC-SINE,  UNMOD):      Cover- 
age Fraction  and  Averase  Half  Lensth. 
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Figure  1: 


CONFIDENCE  INTERVAL(80%,  ARC-SINE,  UNMOD):      Cover- 
age Fraction  and  Average  Half  Length. 
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Figure  16. 


CONFIDENCE  INTERVAL(90%,  ARC-SINE,  MOD  1): 

Fraction  and  Average  Half  Length. 
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Figure   17.       CONFIDENCE  INTERVAL(80%,  ARC-SINE,  MOD  1):      Coverage 


Fraction  and  Average  Half  Length. 
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IV.     CONCLUSIONS 

This  thesis  considers  the  problem  of  estimating  the  survival  probability  P{D  >  i]  for 
the  first  passage  time  to  state  0  for  a  semi-Markov  process  using  censored  data.  Simu- 
lation is  used  to  study  the  small  sample  behavior  of  three  estimators  and  their  confidence 
interval  procedures. 

One  of  the  estimators  studied  is  the  Kaplan- Meier  estimator  of  the  first  passage 
times  to  state  0.  Both  the  unmodified  Kaplan-Meier  estimator  and  two  modifications, 
MOD  1  and  MOD  2  making  the  estimated  distribution  honest  are  considered.  Another 
estimator  is  the  maximum  likelihood  estimator.  A  third  estimator,  the  asymptotic  re- 
newal estimator,  uses  an  exponential  approximation  to  the  survival  function. 

The  following  conclusions  are  drawn  from  the  simulation  experiment. 

1.  The  modified  Kaplan-Meier  estimators  MOD  1  and  MOD  2  using  the  first  passage 
times  to  state  0  have  a  smaller  bias  for  large  times  than  the  unmodified  K.M.E.. 
However,  in  the  medium  range  of  times,  the  two  modified  procedures  MOD  1  and 
MOD  2  have  a  larger  bias  than  that  of  the  unmodified  K.M.E..  MOD  1  is  slightly 
better  than  MOD  2. 

2.  Modifying  the  Kaplan-Meier  estimates  of  the  sojourn  time  distributions  in  the 
asymptotic  renewal  estimate  does  not  improve  its  performance. 

3.  The  asymptotic  normal  confidence  intervals  for  the  Kaplan-Meier  estimator  of  the 
first  passage  times  to  state  0  using  the  arc-sine  transformation  have  a  slightly  better 
coverage  than  those  using  the  log  transformation. 

4.  The  confidence  intervals  using  the  log  transformed  estimators  are  preferred  to 
those  using  the  arc-sine  transformed  estimators  for  the  maximum  likelihood  esti- 
mator and  the  asymptotic  renewal  estimator. 
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5.  The  confidence  intervals  for  the  jackknifed  log  transformed  asymptotic  renewal 
estimator  have  the  correct  coverage  for  all  but  the  smallest  times.  This  estimator 
makes  no  assumptions  concerning  the  parametric  form  of  the  sojourn  time  distrib- 
utions. It's  expected  that  it  will  also  perform  well  in  cases  in  which  the  sojourn  time 
distributions  are  not  exponential. 

6.  The  confidence  intervals  for  the  maximum  likelihood  estimator  have  the  correct 
coverage,  and  also  have  the  smallest  average  half  length.  However  the  estimators 
depend  on  the  parametric  form  of  the  sojourn  time  distributions.  If  the  parametric 
form  is  incorrectly  specified,  the  M.L.E.  can  be  quite  biased,  Kim[Ref.  1]. 
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APPENDIX  A.      TRUE  SURVIVAL  PROBABILITY 


Table  1.     TRUE  SURVIVAL  PROBABILITY  FOR  MODEL 


TIME 

P{D  >  t) 

1.00 

0.6644 

3.00 

0.3554 

5.00 

0.1974 

7.00 

0.1099 

9.00 

0.0612 

11.00 

0.0340 

13.00 

0.0189 

15.00 

0.0105 
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APPENDIX  B.      NUMBER  OF  K.M.E.  DEFINED  AT  TIME  T 


Table  2.     NUMBER  OF  KAPLAN-MEIER  ESTIMATES  DEFINED  AT  TIME  T: 

C=l/2 


TIME 

KMO 

KM  I 

KM2 

1.00 

500 

500 

496 

3.00 

496 

428 

369 

5.00 

366 

314 

334 

7.00 

215 

308 

333 

9.00 

171 

308 

333 

11.00 

165 

308 

333 

13.00 

164 

308 

333 

15.00 

164 

308 

333 

Table  3.     NUMBER  OF  KAPLAN-MEIER  ESTIMATES  DEFINED  AT  TIME  T: 
C=  1/5 


TIME 

KMO 

KM1 

KM2 

1.00 

500 

500 

500 

3.00 

500 

488 

463 

5.00 

496 

436 

425 

7.00 

456 

421 

419 

9.00 

375 

420 

419 

11.00 

330 

420 

419 

13.00 

309 

420 

419 

15.00 

305 

420 

419 
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Table  4.     NUMBER  OF  KAPLAN-MEIER  ESTIMATES  DEFINED  AT  TIME  T: 

C=l/10 


TIME 

KMO 

KM1 

KM2 

1.00 

500 

500 

500 

3.00 

500 

497 

484 

5.00 

499 

462 

455 

7.00 

489 

451 

451 

9.00 

466 

447 

450 

11.00 

436 

447 

450 

13.00 

412 

447 

450 

15.00 

399 

447 

450 

Table  5.     NUMBER  OF  KAPLAN-MEIER  ESTIMATES  DEFINED  AT  TIME  T: 

C=  1/100 


TIME 

KMO 

KM1 

KM2 

1.00 

500 

500 

500 

3.00 

500 

500 

499 

5.00 

500 

498 

494 

7.00 

500 

494 

493 

9.00 

499 

494 

493 

11.00 

499 

494 

493 

13.00 

498 

494 

493 

15.00 

494 

494 

493 
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APPENDIX  C.      STATISTICS  OF  THREE  ESTIMATORS 


Table  6.     STATISTICS  OF  THREE  ESTIMATORS  :  UNMOD,  C=  1/2 


TIME 

K.M.E 

M.L.E 

A.R.E 

MEAN 

SRMSE 

MEAN 

SRMSE 

MEAN 

SRMSE 

1.0 

3.0 

5.0 

7.0 

9.0 

11.0 

13.0 

15.0 

0.6608 
0.3672 
0.2200 
0.1778 
0.1700 
0.1683 
0.1683 
0.1683 

0.0057 
0.0124 
0.0198 
0.0263 
0.0342 
0.0403 
0.0446 
0.0472 

0.6654 
0.3633 
0.2098 
0.1237 

0.0742 
0.0452 
0.0279 
0.0175 

0.0033 
0.0064 
0.0055 
0.0037 
0.0023 
0.0013 
0.0007 
0.0004 

0.5758 
0.2739 
0.1370 
0.0713 
0.0383 
0.0212 
0.0120 
0.0070 

0.0196 
0.0178 
0.0099 
0.0045 
0.0019 
0.0008 
0.0003 
0.0002 

Table  7.     STATISTICS  OF  THREE  ESTIMATORS  :  UNMOD,  C=  1/5 


TIME 

K.M.E 

M.L.E 

A.R.E 

MEAN 

SRMSE 

MEAN 

SRMSE 

MEAN 

SRMSE 

1.0 

3.0 

5.0 

7.0 

9.0 

11.0 

13.0 

15.0 

0.6688 
0.3540 
0.1982 
0.1135 
0.0708 
0.0507 
0.0435 
0.0399 

0.0048 
0.0068 
0.0063 
0.0054 
0.0050 
0.0043 
0.0042 
0.0042 

0.6599 
0.3547 
0.2004 
0.1148 
0.0665 
0.0388 
0.0229 
0.0136 

0.0024 
0.0042 
0.0033 
0.0020 
0.0010 
0.0005 
0.0003 
0.0001 

0.6373 
0.3389 
0.1838 
0.1014 
0.0568 
0.0323 
0.0185 
0.0 108 

0.0043 
0.0056 
0.0040 
0.0022 
0.0011 
0.0005 
0.0002 
0.0001 
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Table  8.     STATISTICS  OF  THREE  ESTIMATORS  :  UNMOD,  C=  1/10 


TIME 

K.M.E 

M.L.E 

A.R.E 

MEAN 

SRMSE 

MEAN 

SRMSE 

MEAN 

SRMSE 

1.0 

3.0 

5.0 

7.0 

9.0 

11.0 

13.0 

15.0 

0.6628 
0.3518 
0.1962 
0.1087 
0.0644 
0.0380 
0.0249 
0.0179 

0.0046 
0.0050 
0.0041 
0.0032 
0.0023 
0.0017 
0.0013 
0.0011 

0.6604 
0.3535 
0.1981 
0.1124 
0.0643 
0.0370 
0.0215 
0.0126 

0.0021 
0.0035 
0.0026 
0.0015 
0.0008 
0.0004 
0.0002 
0.0001 

0.6428 
0.3476 
0.1905 
0.1057 
0.0593 
0.0336 
0.0192 
0.0111 

0.0031 
0.0040 
0.0029 
0.0016 
0.0008 
0.0004 
0.0002 
0.0001 

Table  9.     STATISTICS  OF  THREE  ESTIMATORS  :  UNMOD,  C=  1/100 


TIME 

K.M.E 

M.L.E 

A.R.E 

MEAN 

SRMSE 

MEAN 

SRMSE 

MEAN 

SRMSE 

1.0 

3.0 

5.0 

7.0 

9.0 

11.0 

13.0 

15.0 

0.6633 
0.3558 
0.1974 
0.1078 
0.0593 
0.0327 
0.0179 
0.0110 

0.0044 
0.0045 
0.0034 
0.0019 
0.0011 
0.0007 
0.0004 
0.0003 

0.6591 
0.3513 
0.1956 
0.1100 
0.0623 
0.0355 
0.0204 
0.0118 

0.0015 
0.0027 
0.0021 
0.0012 
0.0006 
0.0003 
0.0001 
0.0001 

0.6434 
0.3530 
0.1954 
0.1090 
0.0613 
0.0347 
0.0198 
0.0114 

0.0025 
0.0027 
0.0021 
0.0012 
0.0006 
0.0003 
0.0001 
0.0001 
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Table  10.     STATISTICS  OF  THREE  ESTIMATORS  :  MOD  1,  C=  1/2 


TIME 

K.M.E 

M.L.E 

A.R.E 

MEAN 

SRMSE 

MEAN 

SRMSE 

MEAN 

SRMSE 

1.0 

3.0 

5.0 

7.0 

9.0 

11.0 

13.0 

15.0 

0.6608 
0.3659 
0.1521 
0.0355 
0.0050 
0.0002 
0.0000 
0.0000 

0.0057 
0.0131 
0.0259 
0.0149 
0.0047 
0.0012 
0.0004 
0.0001 

0.6654 
0.3633 
0.2098 
0.1237 
0.0742 
0.0452 
0.0279 
0.0175 

0.0033 
0.0064 
0.0055 
0.0037 
0.0023 
0.0013 
0.0007 
0.0004 

0.4914 
0.2419 
0.1239 
0.0656 
0.0358 
0.0200 
0.0115 
0.0067 

0.0665 
0.0293 
0.0127 
0.0052 
0.0021 
0.0008 
0.0003 
0.0002 

Table  11.     STATISTICS  OF  THREE  ESTIMATORS  :  MOD  I,  C=  1/5 


TIME 

K.M.E 

M.L.E 

A.R.E 

MEAN 

SRMSE 

MEAN 

SRMSE 

MEAN 

SRMSE 

1.0 

3.0 

5.0 

7.0 

9.0 

11.0 

13.0 

15.0 

0.6688 
0.3540 
0.1973 
0.1046 
0.0465 
0.0174 
0.0055 
0.0011 

0.0048 
0.0068 
0.0065 
0.0063 
0.0051 
0.0024 
0.0009 
0.0002 

0.6599 
0.3547 
0.2004 
0.1148 
0.0665 
0.0388 
0.0229 
0.0136 

0.0024 
0.0042 
0.0033 
0.0020 
0.0010 
0.0005 
0.0003 
0.0001 

0.6068 
0.3247 
0.1770 
0.0981 
0.0551 
0.0314 
0.0181 
0.0105 

0.0142 
0.0086 
0.0049 
0.0024 
0.0011 
0.0005 
0.0002 
0.0001 
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Table   12.     STATISTICS  OF  THREE  ESTIMATORS  :  MOD  1,  C=  1/10 


TIME 

K.M.E 

M.L.E 

A.R.E 

MEAN 

SRMSE 

MEAN 

SRMSE 

MEAN 

SRMSE 

1.0 

0.6628 

0.0046 

0.6604 

0.0021 

0.6247 

0.0086 

3.0 

0.3518 

0.0050 

0.3535 

0.0035 

0.3389 

0.0056 

5.0 

0.1962 

0.0041 

0.1981 

0.0026 

0.1862 

0.0034 

7.0 

0.1074 

0.0034 

0.1124 

0.0015 

0.1035 

0.0018 

9.0 

0.0605 

0.0025 

0.0643 

0.0008 

0.0582 

0.0008 

11.0 

0.0304 

0.0016 

0.0370 

0.0004 

0.0330 

0.0004 

13.0 

0.0149 

0.0009 

0.0215 

0.0002 

0.0189 

0.0002 

15.0 

0.0060 

0.0004 

0.0126 

0.0001 

0.0109 

0.0001 

Table  13.     STATISTICS  OF  THREE  ESTIMATORS  :  MOD  1,  C=  1/100 


TIME 

K.M.E 

M.L.E 

A.R.E 

MEAN 

SRMSE 

MEAN 

SRMSE 

MEAN 

SRMSE 

1.0 

3.0 

5.0 

7.0 

9.0 

11.0 

13.0 

15.0 

0.6633 
0.3558 
0.1974 
0.1078 
0.0592 
0.0326 
0.0178 
0.0108 

0.0044 
0.0045 
0.0034 
0.0019 
0.0011 
0.0007 
0.0004 
0.0003 

0.6591 
0.3513 
0.1956 
0.1100 
0.0623 
0.0355 
0.0204 
0.0118 

0.0015 
0.0027 
0.0021 
0.0012 
0.0006 
0.0003 
0.0001 
0.0001 

0.6419 
0.3522 
0.1950 
0.1088 
0.0612 
0.0347 
0.0198 
0.0114 

0.0027 
0.0028 
0.0021 
0.0012 
0.0006 
0.0003 
0.0001 
0.0001 
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Table   14.     STATISTICS  OF  THREE  ESTIMATORS  :  MOD  2,  C=  1/2 


TIME 

K.M.E 

M.L.E 

A.R.E 

MEAN 

SRMSE 

MEAN 

SRMSE 

MEAN 

SRMSE 

1.0 

0.6608 

0.0057 

0.6654 

0.0033 

0.4966 

0.0621 

3.0 

0.3426 

0.0162 

0.3633 

0.0064 

0.2439 

0.0284 

5.0 

0.1050 

0.0233 

0.2098 

0.0055 

0.1247 

0.0125 

7.0 

0.0205 

0.0116 

0.1237 

0.0037 

0.0659 

0.0052 

9.0 

0.0032 

0.0041 

0.0742 

0.0023 

0.0359 

0.0021 

11.0 

0.0002 

0.0012 

0.0452 

0.0013 

0.0201 

0.0008 

13.0 

0.0000 

0.0004 

0.0279 

0.0007 

0.0115 

0.0003 

15.0 

0.0000 

0.0001 

0.0175 

0.0004 

0.0067 

0.0002 

Table  15.     STATISTICS  OF  THREE  ESTIMATORS  :  MOD  2,  C=  1/5 


TIME 

K.M.E 

M.L.E 

A.R.E 

MEAN 

SRMSE 

MEAN 

SRMSE 

MEAN 

SRMSE 

1.0 

3.0 

5.0 

7.0 

9.0 

11.0 

13.0 

15.0 

0.6688 
0.3540 
0.1946 
0.0984 
0.0420 
0.0161 
0.0054 
0.0009 

0.0048 
0.0068 
~  0.0068 
0.0061 
0.0046 
0.0022 
0.0009 
0.0002 

0.6599 
0.3547 
0.2004 
0.1148 
0.0665 
0.0388 
0.0229 
0.0136 

0.0024 
0.0042 
0.0033 
0.0020 
0.0010 
0.0005 
0.0003 
0.0001 

0.6079 
0.3252 
0.1772 
0.0982 
0.0552 
0.0314 
0.0181 
0.0105 

0.0135 
0.0084 
0.0049 
0.0024 
0.0011 
0.0005 
0.0002 
0.0001 
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Table   16.     STATISTICS  OF  THREE  ESTIMATORS  :  MOD  2,  C=  1/10 


TIME 

K.M.E 

M.L.E 

A.R.E 

MEAN 

SRMSE 

MEAN 

SRMSE 

MEAN 

SRMSE 

1.0 

3.0 

5.0 

7.0 

9.0 

11.0 

13.0 

15.0 

0.6628 
0.3518 
0.1959 
0.1067 
0.0594 
0.0294 
0.0140 
0.0057 

0.0046 
0.0050 
0.0041 
0.0034 
0.0024 
0.0015 
0.0008 
0.0003 

0.6604 
0.3535 
0.1981 
0.1124 
0.0643 
0.0370 
0.0215 
0.0126 

0.0021 
0.0035 
0.0026 
0.0015 
0.0008 
0.0004 
0.0002 
0.0001 

0.6254 
0.3392 
0.1864 
0.1036 
0.0582 
0.0330 
0.0189 
0.0109 

0.0081 
0.0055 
0.0034 
0.0017 
0.0008 
0.0004 
0.0002 
0.0001 

Table   17.     STATISTICS  OF  THREE  ESTIMATORS  :  MOD  2,  C=  1/100 


TIME 

K.M.E 

M.L.E 

A.R.E 

MEAN 

SRMSE 

MEAN 

SRMSE 

MEAN 

SRMSE 

1.0 

0.6633 

0.0044 

0.6591 

0.0015 

0.6419 

0.0027 

3.0 

0.3558 

0.0045 

0.3513 

0.0027 

0.3522 

0.0028 

5.0 

0.1974 

0.0034 

0.1956 

0.0021 

0.1950 

0.0021 

7.0 

0.1078 

0.0019 

0.1100 

0.0012 

0.1088 

0.0012 

9.0 

0.0592 

0.0011 

0.0623 

0.0006 

0.0612 

0.0006 

11.0 

0.0326 

0.0007 

0.0355 

0.0003 

0.0347 

0.0003 

13.0 

0.0178 

0.0004 

0.0204 

0.0001 

0.0198 

0.0001 

15.0 

0.0108 

0.0003 

0.0118 

0.0001 

0.0114 

0.0001 
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APPENDIX  D.     CONFIDENCE  INTERVALS 


Table  18.     TWO-SIDED  90  ( 

Vo  COVERAGE  FRACTION(UNMOD,  C  = 

1/2) 

TIME 

COV- 
ERAGE 

K.M.E 

M.L.E 

A.R.E 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

1.0 

toohigh 

44(.09) 

30(.06) 

42(.08) 

26(.05) 

12(.02) 

11(.02) 

cover 

444(.89) 

453(.91) 

445(.89) 

455(.91) 

434(.87) 

401(.80) 

too  low 

12(.02) 

17(.03) 

13(.03) 

19(.04) 

54(.ll) 

88(.18) 

3.0 

too  high 

55(.ll) 

39(.08) 

35(.07) 

29(.06) 

12(.02) 

10(.02) 

cover 

439(.8S) 

43S(.S8) 

449(.90) 

447(.89) 

453(.91) 

402(.80) 

too  low 

6(.01) 

23(.05) 

16(.03) 

24(.05) 

35(.07) 

88(.18) 

5.0 

too  high 

100(.20) 

78(.16) 

39(.08) 

27(.05) 

12(.02) 

9(.02) 

cover 

321(.64) 

340(.68) 

444(,89) 

448(.90) 

452(.90) 

382(.76) 

too  low 

79(.16) 

82(.16) 

17(.03) 

25(.05) 

36(.07) 

109(.22) 

7.0 

too  high 

162(.32) 

134(.27) 

41(.08) 

26(.05) 

12(.02) 

8(.02) 

cover 

194(39) 

222(.44) 

443(.89) 

448(.90) 

453(.91) 

355(.71) 

too  low 

144(.29) 

144(.29) 

16(.03) 

26(.05) 

35(.07) 

137(.27) 

9.0 

too  high 

224(.45) 

190(.38) 

43(.09) 

25(.05) 

12(.02) 

8(.02) 

cover 

115(.23) 

149(.30) 

440(.88) 

440(.88) 

452(.90) 

344(.69) 

too  low 

161(.32) 

161(.32) 

17(.03) 

35(.07) 

36(.07) 

148(.30) 

11.0 

too  high 

2S5(.57) 

217(.43) 

46(.09) 

24(.05) 

12(.02) 

8(.02) 

cover 

51(.10) 

119(.24) 

437(.S7) 

435(.87) 

454(.91) 

322{.64) 

too  low 

164(.33) 

164(.33) 

17(.03) 

41(.0S) 

34(.07) 

170(.34) 

13.0 

too  high 

323(.65) 

250(.50) 

46(.09) 

19(.04) 

12(.02) 

9(.02) 

cover 

13(.03) 

86(.17) 

437(.87) 

437(.87) 

454(.91) 

305(.61) 

too  low 

164(.33) 

164(.33) 

17(.03) 

44(.09) 

34(.07) 

186(.37) 

15.0 

too  high 

336(.67) 

290{.58) 

47(.09) 

16(.03) 

12(.02) 

9(.02) 

cover 

0(.00) 

46(.09) 

436(.87) 

434(.87) 

454(.91) 

2S5(.57) 

too  low 

164(.33) 

164(.33) 

17(.03) 

50(.10) 

34(.07) 

206(.41) 
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Table  19.     AVERAGE  AND  STANDARD  DEVIATION  OF  HALF  LENGTH(90% 
CI,  UNMOD,  C=l/2) 


TIME 

K.M.E 

M.L.E 

A.R.E 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

1.0 

.1245   .1223 

.0959   .0949 

.2032  .1862 

(.0107)  (.0102) 

(.0080)  (.0077) 

(.0841)  (.0719) 

3.0 

.1783   .1647 

.1337   .1289 

.2583   .1948 

(.0341)  (.0249) 

(.0108)  (.0108) 

(.1215)  (.0824) 

5.0 

.1775   .1494 

.1251   .1164 

.2572   .1490 

(.0904)  (.0722) 

(.0223)  (.0223) 

(.1506)  (.0840) 

7.0 

.1535   .1267 

.1045   .0927 

.2415   .1031 

(.1083)  (.0864) 

(.0295)  (.0287) 

(.1704)  (.0778) 

9.0 

.1463   .1207 

.0830  .0695 

.2243   .0683 

(.1108)  (.0887) 

(.0320)  (.0299) 

(.1831)  (.0660) 

11.0 

.1452   .1197 

.0644   .0505 

.2091   .0444 

(.1113)  (.0891) 

(.0316)  (.0281) 

(.1914)  (.0531) 

13.0 

.1452   .1197 

.0495   .0361 

.1966   .0285 

(.1113)  (.0891) 

(.0295)  (.0249) 

(.1969)  (.0414) 

15.0 

.1452   .1197 

.0379   .0256 

.1866   .0183 

(.1113)  (.0891) 

(.0267)  (.0213) 

(.2009)  (.0316) 

49 


Table  20.     TWO-SIDED  80  %  COVERAGE  FRACTION(UNMOD,  C=  1/2) 


TIME 

COV- 
ERAGE 

K.M.E 

M.L.E 

A.R.E 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

1.0 

too  high 

68(.14) 

59(.12) 

61(.12) 

54(.ll) 

33(.07) 

29(.06) 

cover 

406(.81) 

396(.79) 

407(.81) 

405(.81) 

374(.75) 

342(.68) 

too  low 

26(.05) 

45(.09) 

32(.06) 

41(.08) 

93(.19) 

129(.26) 

3.0 

too  high 

94(.19) 

75(.15) 

62(.12) 

47(.09) 

29(.06) 

26(.05) 

cover 

381(.76) 

382(.76) 

404(.80) 

409(.82) 

401(.80) 

343(.69) 

too  low 

25(.05) 

43(.09) 

34(.07) 

44(.09) 

70(.14) 

131(.26) 

5.0 

too  high 

129(.26) 

108(.22) 

64(.13) 

48(.10) 

25(.05) 

19(.04) 

cover 

291(.58) 

302(.60) 

406(.81) 

405(.81) 

407(.81) 

328(.76) 

too  low 

80(.16) 

90(.1S) 

30(.06) 

47(.09) 

68(.14) 

153(.31) 

7.0 

too  high 

191(.3S) 

164(.33) 

65(.13) 

48(.10) 

24(.05) 

17(.03) 

cover 

165(.33) 

192(.38) 

406(.81) 

398(.80) 

403(.81) 

307(.61) 

too  low 

144(.29) 

144(.29) 

29(.06) 

54(.ll) 

73(.15) 

176(.35) 

9.0 

too  high 

257(.51) 

217{.43) 

67(.13) 

47(.09) 

22(.04) 

16(.03) 

cover 

82(.16) 

122(.24) 

405(.81) 

396(.79) 

403(.81) 

289(.58) 

too  low 

161(.32) 

161(.32) 

28(.06) 

57(.ll) 

75(.15) 

195(.39) 

11.0 

too  high 

307(.61) 

267(.53) 

66(.13) 

45(.09) 

21(.04) 

15(.03) 

cover 

29(.06) 

69(.14) 

405(.81) 

392(.78) 

403(.81) 

273(.55) 

too  low 

164(.33) 

164(.33) 

29(.06) 

63(.13) 

76(.15) 

212(.42) 

13.0 

too  high 

334(.67) 

308(.62) 

69(.14) 

39(.08) 

21  (.04) 

12(.02) 

cover 

2(.00) 

28(.06) 

402(.80) 

396(.79) 

404(.81) 

267(.53) 

too  low 

164(.33) 

164(.33) 

29(.06) 

65(.13) 

75(.15) 

221(.44) 

15.0 

too  high 

336(.67) 

333(.67) 

70(.14) 

34(.07) 

21(.04) 

12(.02) 

cover 

0(.00) 

3(.01) 

400(.80) 

400(.80) 

404(.81) 

252(.50) 

too  low 

164(.33) 

164(.33) 

30(.06) 

66(.13) 

75(.15) 

236(.47) 

50 


Table  21.     AVERAGE  AND  STANDARD  DEVIATION  OF  HALF  LENGTH(80% 
C.I,  UNMOD,  C=l/2) 


TIME 

K.M.E 

M.L.E 

A.R.E 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

1.0 

.0968   .0957 

.0746  .0742 

.1596   .1447 

(.0082)  (.0080) 

(.0062)  (.0061) 

(.0748)  (.0602) 

3.0 

.1359   .1296 

.1032   .1010 

.2048   .1510 

(.0240)  (.0200) 

(.0085)  (.0085) 

(.1153)  (.0677) 

5.0 

.1311   .1181 

.0953   .0912 

.2010   .1152 

(.0655)  (.0574) 

(.0175)  (.0175) 

(.14300  (.0677) 

7.0 

.1126   .1002 

.0781   .0726 

.1845   .0795 

(.0784)  (.0686) 

(.0229)  (.0226) 

(.1598)  (.0621) 

9.0 

.1073   .0955 

.0607   .0545 

.1677   .0526 

(.0804)  (.0704) 

(.0245)  (.0236) 

(.1701)  (.0522) 

11.0 

.1065   .0947 

.0459   .0396 

.1534   .0340 

(.0808)  (.0707) 

(.0237)  (.0221) 

(.1765)  (.0417) 

13.0 

.1065   .0947 

.0343   .0283 

.1420   .0218 

(.0808)  (.0707) 

(.0217)  (.0196) 

(.1805)  (.0322) 

15.0 

.1065   .0947 

.0256   .0200 

.1330   .0140 

(.0808)  (.0707) 

(.0192)  (.0167) 

(.1832)  (.0245) 

51 


Table  22.     TWO-SIDED  90  %  COVERAGE  FRACTION(UNMOD,  C=  1/10) 


TIME 

COV- 
ERAGE 

K.M.E 

M.L.E 

A.R.E 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

1.0 

too  high 

34(.07) 

26(.05) 

31(.06) 

24(.05) 

16(.03) 

12(.02) 

cover 

453(.91) 

454(.91) 

448(.90) 

443(.89) 

436(.87) 

428(.86) 

too  low 

13(.03) 

20(.04) 

21  (.04) 

33(.07) 

48(.10) 

60(.12) 

3.0 

too  high 

24(.05) 

16(.03) 

27(.05) 

21(.04) 

24(.05) 

15(.03) 

cover 

460(.92) 

456(.91) 

450(.90) 

444(.89) 

460(.92) 

450(.90) 

too  low 

16(.03) 

28(.06) 

23(.05) 

35(.07) 

16(.03) 

35(.07) 

5.0 

too  high 

37(.07) 

20(.04) 

27(.05) 

19(.04) 

23(.05) 

13(.03) 

cover 

453(.91) 

452(.90) 

452(.90) 

445(.89) 

457(.91) 

451(.90) 

too  low 

10(.02) 

28(.06) 

21(.04) 

36(.07) 

20(.04) 

36(.07) 

7.0 

too  high 

37(.07) 

25(.05) 

27(.05) 

18(.04) 

23(.05) 

12(.02) 

cover 

43S(.8S) 

434(.87) 

454(.91) 

442(.88) 

456(.91) 

444(.89) 

too  low 

25(.05) 

41(.08) 

19(.04) 

40(.08) 

21(.04) 

44(.09) 

9.0 

too  high 

47(.09) 

27(.05) 

26(.05) 

16(.03) 

24(.05) 

9(.02) 

cover 

357(.71) 

377(.75) 

455(.91) 

438(.88) 

455(.91) 

438(.88) 

too  low 

96(.19) 

96(.19) 

19(.04) 

46(.09) 

21(.04) 

53(.ll) 

11.0 

too  high 

59(.12) 

29(.06) 

26(.05) 

12(.02) 

24(.05) 

7(.01) 

cover 

233(.47) 

263(.53) 

455(.91) 

434(.87) 

455(.91) 

428(.86) 

too  low 

208(.42) 

208(.42) 

19(.04) 

54C.11) 

21(.04) 

65(.13) 

13.0 

too  high 

59(.12) 

38(.08) 

26(.05) 

11  (.02) 

24(.05) 

6(.01) 

cover 

157(.31) 

175(.35) 

454(.91) 

433(.87) 

455(.91) 

422(.84) 

too  low 

287(.57) 

287(.57) 

20(.04) 

56(.ll) 

21(.04) 

72(.14) 

15.0 

too  high 

83(.17) 

35(.07) 

26(.05) 

7(.01) 

24(.05) 

4(.01) 

cover 

73(.15) 

121(.24) 

456(.91) 

436(.87) 

455(.91) 

419(.84) 

too  low 

344(.69) 

344(.69) 

18(.04) 

57(.ll) 

21  (.04) 

77(.15) 

52 


Table  23.     AVERAGE  AND  STANDARD  DEVIATION  OF  HALF  LENGTH(90% 
C.I,  UNMOD,  C=l/10) 


TIME 

K.M.E 

M.L.E 

A.R.E 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

1.0 

.1121   .1105 

.0718   .0714 

.0932   .0903 

(.0064)  (.0062) 

(.0061)  (.0060) 

(.0317)  (.0279) 

3.0 

.1223   .1185 

.0972   .0953 

.1212  .1100 

(.0062)  (.0066) 

(.0042)  (.0044) 

(.0421)  (.0297) 

5.0 

.1123   .1046 

.0871   .0837 

.1160   .0969 

(.0132)  (.0142) 

(.0101)  (.0103) 

(.0522)  (.0299) 

7.0 

.0988   .0847 

.0688   .0644 

.0981   .0738 

(.0285)  (.0267) 

(.0134)  (.0135) 

(.0576)  (.0280) 

9.0 

.0806   .0640 

.0512   .0464 

.0789   .0523 

(.0437)  (.0363) 

(.0139)  (.0137) 

(.0596)  (.0241) 

11.0 

.0573   .0432 

.0369   .0321 

.0624   .0355 

(.0513)  (.0398) 

(.0128)  (.0122) 

(.0601)  (.0195) 

13.0 

.0416   .0305 

.0261   .0217 

.0494   .0235 

(.0506)  (.0380) 

(.0109)  (.0101) 

(.0601)  (.0151) 

15.0 

.0306   .0223 

.0182   .0145 

.0395   .0153 

(.0472)  (.0350) 

(.0089)  (.0079) 

(.0602)  (.0113) 

53 


Table  24.     TWO-SIDED  80  %  COVERAGE  FRACTION(UNMOD,  C=  1/10) 


TIME 

COV- 
ERAGE 

K.M.E 

M.L.E 

A.R.E 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

1.0 

too  high 

61(.12) 

47(.09) 

56(.ll) 

48(.10) 

38(.08) 

31(.06) 

cover 

404(.81) 

402(.80) 

400(.80) 

403(.81) 

382(.76) 

379(.76) 

too  low 

35(.07) 

51(10) 

44(.09) 

49(.10) 

80(.16) 

90(.18) 

3.0 

too  high 

49(.10) 

38(.08) 

52(.10) 

42(.08) 

57(.ll) 

43(.09) 

cover 

415(.83) 

417(.83) 

401(.80) 

403(.81) 

401(.80) 

398(.80) 

too  low 

36(.07) 

45(.09) 

47(.09) 

55(.ll) 

42(.08) 

59(.12) 

5.0 

too  high 

56(.ll) 

47(.09) 

52(.10) 

40(.08) 

53(.ll) 

40(.08) 

cover 

41 1(.S2) 

392(.78) 

397(.79) 

395(.79) 

406(.81) 

395(.79) 

too  low 

33(.07) 

61(.12) 

51(.10) 

65(.13) 

41(.08) 

65(.13) 

7.0 

too  high 

55(.ll) 

43(.09) 

51(.10) 

39(.08) 

50(.01) 

35(.07) 

cover 

410(.S2) 

390{.78) 

395(.79) 

396(.79) 

407(.81) 

392(.78) 

too  low 

35(.07) 

67(.13 

54(.ll) 

65(.13) 

43(.09) 

73(.15) 

9.0 

too  high 

7l(.14) 

53(.ll) 

51(.10) 

34(.07) 

50(.10) 

32(.06) 

cover 

333(.67) 

351(.70) 

394(.79) 

399(.80) 

409(.82) 

387(.77) 

too  low 

96(.19) 

96(.19) 

55(.12) 

67(.13) 

41(.08) 

8 1  ( .  1 6 ) 

11.0 

too  high 

79(.16) 

59(.12) 

53(.ll) 

31  (.06) 

49(.10) 

26(.05) 

cover 

213(.43) 

233(.47) 

392(.78) 

400(.80) 

408(.82) 

382(.76) 

too  low 

208(.42) 

208(.42) 

55(.ll) 

69(.14) 

43(.09) 

92(.18) 

13.0 

too  high 

84(.17) 

56(1 1) 

55(.ll) 

29(.06) 

49(.10) 

25(.05) 

cover 

129(.26) 

157(.31) 

391(.78) 

398(.89) 

406(.81) 

379(.76) 

too  low 

287(.57) 

287(.57) 

54(.ll) 

73(.15) 

45(.09) 

96(.19) 

15.0 

too  high 

110(.22) 

58(.12) 

54(.ll) 

26(.05) 

50(.10) 

21(.04) 

cover 

46(.09) 

98(.20) 

392(.78) 

398(.80) 

405(.81) 

375(.75) 

too  low 

344(.69) 

844(.69) 

54(.ll) 

76(.15) 

45(.09) 

104(.21) 

54 


Table  25.     AVERAGE  AND  STANDARD  DEVIATION  OF  HALF  LENGTH(80% 
C.I,  UNMOD,  C=  1/10) 


TIME 

K.M.E 

M.L.E 

A.R.E 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

1.0 

.0872  .0864 

.0559   .0557 

.0701   .0683 

(.0049)  (.0048) 

(.0047)  (.0047) 

(.0238)  (.0213) 

3.0 

.0945   .0927 

.0754   .0745 

.0905   .0834 

(.0050)  (.0051) 

(.0033)  (.0035) 

(.0309)  (.0227) 

5.0 

.0855   .0819 

.0670   .0654 

.0851   .0736 

(.0106)  (.0112) 

(.0080)  (.0081) 

(.0369)  (.0230) 

7.0 

.0729   .0664 

.0524   .0503 

.0704   .0560 

(.0216)  (.0209) 

(.0105)  (.0106) 

(.0407)  (.0214) 

9.0 

.0578   .0502 

.0385   .0362 

.0551   .0397 

(.0317)  (.0285) 

(.0108)  (.0107) 

(.0430)  (.0184) 

11.0 

.0403   .0339 

.0273   .0251 

.0422   .0269 

(.0364)  (.0312) 

(.0098)  (.0095) 

(.0432)  (.0149) 

13.0 

.0290   .0240 

.0190  .0169 

.0322   .0178 

(.0355)  (.0298) 

(.0083)  (.0079) 

(.0430)  (.0115) 

15.0 

.0213  ,  .0175 

.0130   .0113 

.0247   .0116 

(.0330)  (.0275) 

(.0067)  (.0062) 

(.0427)  (.0085) 

55 


Table  26.     TWO-SIDED  90  %  COVERAGE  FRACTION(MOD  1,  C=  1/2) 


TIME 

COV- 
ERAGE 

K.M.E 

M.L.E 

A.R.E 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

1.0 

too  high 

44(.09) 

30(.06) 

42(.08) 

26(.05) 

9(.02) 

8(.02) 

cover 

444(.89) 

453(.91) 

445(.89) 

455(.91) 

451(.90) 

377(.75) 

too  low 

12(.02) 

17(.03) 

13(.03) 

19(.04) 

40(.08) 

115(.24) 

3.0 

too  high 

55(.ll) 

39(.08) 

35(.07) 

29(.06) 

9(.02) 

7(.01) 

cover 

439(.88) 

438(.88) 

449(.90) 

447(.89) 

457(.91) 

382(.76) 

too  low 

6(.01) 

23(.05) 

16(.03) 

24(.05) 

34(.07) 

111(.22) 

5.0 

too  high 

100(.20) 

78(.16) 

39(.08) 

27(.05) 

9(.02) 

10(.02) 

cover 

321(.64) 

340(.68) 

444(.89) 

448(.90) 

460(.92) 

354(.71) 

too  low 

79(.16) 

82(.16) 

17(.03) 

25(.05) 

31(.06) 

136(.27) 

7.0 

too  high 

162(.32) 

134(.27) 

41(.08) 

26(.05) 

10(.02) 

10(.02) 

cover 

194(.39) 

222(.44) 

443(.89) 

448(.90) 

456(.91) 

335(.67) 

too  low 

144(.29) 

144(.29) 

16(.03) 

26(.05) 

34(.07) 

155(.31) 

9.0 

too  high 

224(.45) 

190(38) 

43(.09) 

25(.05) 

10(.02) 

10(.02) 

cover 

115(.23) 

149(.30) 

440(.88) 

440(.8S) 

458(.92) 

324(.65) 

too  low 

161(.32) 

16K.32) 

17(.03) 

35(.07) 

32(.06) 

166(33) 

11.0 

too  high 

285(.57) 

217(.43) 

46(.09) 

24(.05) 

11(.02) 

10(.02) 

cover 

51(.10) 

119(,24) 

437(.87) 

435(.87) 

456(.91) 

307(.61) 

too  low 

164(.33) 

164(.33) 

17(.03) 

41(.08) 

33(.07) 

183(37) 

13.0 

too  high 

323(.65) 

250(.50) 

46(.09) 

19(.04) 

11(.02) 

11(.02) 

cover 

13(.03) 

86(.17) 

437(.87) 

437(.87) 

454(.91) 

290(.58) 

too  low 

164(.33) 

164(.33) 

17(.03) 

44(.09) 

35(.07) 

119(.40) 

15.0 

too  high 

336(.67) 

290(.58) 

47(.09) 

16(.03) 

11(.02) 

12(.02) 

cover 

0(.00) 

46(.09) 

436(.87) 

434(.87) 

454(.91) 

274(.55) 

too  low 

164(.33) 

164(.33) 

17(.03) 

51(.10) 

35(.07) 

214(.43) 

56 


Table  27.     AVERAGE  AND  STANDARD  DEVIATION  OF  HALF  LENGTH(90% 
C.I,  MOD  1,  C=  1/2) 


TIME 

K.M.E 

M.L.E 

A.R.E 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

1.0 

.1245   .1223 

.0959   .0949 

.2871   .2379 

(.0107)  (.0102) 

(.0080)  (.0077) 

(.1289)  (.1119) 

3.0 

.1783   .1647 

.1337   .1289 

.3048   .2088 

(.0341)  (.0249) 

(.0108)  (.0108) 

(.1359)  (.1049) 

5.0 

.1775   .1494 

.1251   .1164 

.2881   .1527 

(.0904)  (.0722) 

(.0223)  (.0223) 

(.1586)  (.1030) 

7.0 

.1535   .1267 

.1045   .0927 

.2657   .1040 

(.1083)  (.0864) 

(.0295)  (.0287) 

(.1775)  (.0902) 

9.0 

.1463   .1207 

.0830   .0695 

.2448   .0684 

(.1108)  (.0887) 

(.0320)  (.0299) 

(.1902)  (.0736) 

9.0 

.1452   .1197 

.0644   .0505 

.2276   .0442 

(.1113)  (.0891) 

(.0316)  (.0281) 

(.1991)  (.0578) 

9.0 

.1452   .1197 

.0495   .0361 

.2137   .0284 

(.1113)  (.0891) 

(.0295)  (.0249) 

(.2050)  (.0444) 

9.0 

.1452   .1197 

.0379   .0256 

.2023   .0182 

(.1113)  (.0891) 

(.0267)  (.0213) 

(.2089)  (.0336) 

57 


Table  28.     TWO-SIDED  80  %  COVERAGE  FRACTION(MOD  1,  C=  1/2) 


TIME 

COV- 
ERAGE 

K.M.E 

M.L.E 

A.R.E 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

1.0 

too  high 

68(.14) 

59(.12) 

61(.12) 

54(.ll) 

28(.06) 

25(.05) 

cover 

406(.81) 

396(.79) 

407(.81) 

405(.81) 

391(.78) 

322(.64) 

too  low 

26(.05) 

45(.09) 

32(.06) 

41(.08) 

81(.16) 

153(.31) 

3.0 

too  high 

94(.19) 

75(.15) 

62(.12) 

47(.09) 

22(.04) 

20(.04) 

cover 

381(.76) 

382(.76) 

404(.81) 

409(.82) 

406(.81) 

318(.64) 

too  low 

25(.05) 

43(.09) 

34(.07) 

44(.09) 

72(.14) 

162(.32) 

5.0 

too  high 

129(.26) 

108(.22) 

64(.13) 

48(.10) 

20(.04) 

18(.04) 

cover 

291(.58) 

302(.60) 

406(.81) 

405(.81) 

409(.82) 

305(.61) 

too  low 

80(.16) 

90(.18) 

30(.06) 

47(.09) 

7K.14) 

177(.35) 

7.0 

too  high 

191(.38) 

164(.33) 

65(.13) 

48(.10) 

21(.04) 

17(.03) 

cover 

1 65(33) 

192(.38) 

406(.81) 

398(.80) 

401(.80) 

292(.58) 

too  low 

144(.29) 

144(.29) 

29(.06) 

54(.ll) 

78(.16) 

191(.38) 

9.0 

too  high 

257(.51) 

217(.43) 

67(.13) 

47(.09) 

19(.04) 

18(.04) 

cover 

82(.16) 

122(.24) 

405(.81) 

396(.79) 

403(.81) 

275(.55) 

too  low 

161(.32) 

161(.32) 

28(.06) 

57(.ll) 

78(.16) 

207(.41) 

11.0 

too  high 

307(.61) 

267(.53) 

66(.13) 

45(.09) 

19(.04) 

17(.03) 

cover 

29(.06) 

69(.14) 

405(.81) 

392(.78) 

403(,81) 

264(.53) 

too  low 

164(.33) 

164(.33) 

29(.06) 

63(.13) 

78(.16) 

219(.44) 

13.0 

too  high 

334(.67) 

308(.62) 

69(.14) 

39(.08) 

19(.04) 

13(.03) 

cover 

2(.00) 

28(.06) 

402(.80) 

396(.79) 

401(.80) 

253(.51) 

too  low 

164(.33) 

164(.33) 

29(.06) 

65(.13) 

80(.16) 

234(.47) 

15.0 

too  high 

336(.67) 

333(.67) 

70(.14) 

34(.07) 

19(.04) 

12(.02) 

cover 

0(.00) 

3  (.01) 

400(.80) 

400(.80) 

403(.81) 

245(.49) 

too  low 

164(.33) 

164(.33) 

30(.06) 

66(.13) 

78(.16) 

243(.49) 
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Table  19.     AVERAGE  AND  STANDARD  DEVIATION  OF  HALF  LENGTH(80% 
C.I,  MOD  1,  C=l/2) 


TIME 

K.M.E 

M.L.E 

A.R.E 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

1.0 

.0968   .0957 

.0746  .0742 

.2321   .1925 

(.0082)  (.0080) 

(.0062)  (.0061) 

(.1207)  (.0972) 

3.0 

.1359   .1296 

.1032   .1010 

.2440   .1648 

(.0240)  (.0200) 

(.0085)  (.0085) 

(.1349)  (.0879) 

5.0 

.1311   .1181 

.0953   .0912 

.2284   .1193 

(.0655)  (.0574) 

(.0175)  (.0175) 

(.1592)  (.0843) 

7.0 

.1126  .1002 

.0781   .0726 

.2069   .0807 

(.0784)  (.0686) 

(.0229)  (.0226) 

(.1748)  (.0726) 

9.0 

.1073   .0955 

.0607   .0545 

.1865   .0528 

(.0804)  (.0704) 

(.0245)  (.0236) 

(.1829)  (.0586) 

11.0 

.1065   .0947 

.0459   .0396 

.1697   .0340 

(.0808)  (.0707) 

(.0237)  (.0221) 

(.1876)  (.0456) 

13.0 

.1065   .0947 

.0343   .0283 

.1566   .0218 

(.0808)  (.0707) 

(.0217)  (.0196) 

(.1906)  (.0347) 

15.0 

.1065   .0947 

.0256  .0200 

.1463   .0139 

(.0808)  (.0707) 

(.0192)  (.0167) 

(.1924)  (.0261) 
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Table  30.     TWO-SIDED  90  %  COVERAGE  FRACTION(MOD  1,  C=  1/10) 


TIME 

COV- 
ERAGE 

K.M.E 

M.L.E 

A.R.E 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

1.0 

too  high 

34{.07) 

26(.05) 

31(.06) 

24(.05) 

14(.03) 

10(.02) 

cover 

453(.91) 

454(.91) 

448(.90) 

443(.89) 

454(.91) 

440(.88) 

too  low 

13(.03) 

20(.04) 

21(.04) 

33(.07) 

32(.06) 

50(.10) 

3.0 

too  high 

24Y.05) 

16(.03) 

27(.05) 

21  (.04) 

20(.04) 

13(.03) 

cover 

460(.92) 

456(.91) 

450(.90) 

444(.89) 

463(.93) 

448(.90) 

too  low 

16(.03) 

28(.06) 

23(.05) 

35(.67) 

17(.03) 

39(.08) 

5.0 

too  high 

37(.07) 

20(.04) 

27(.05) 

19(.04) 

19(.04) 

11(.02) 

cover 

453(.91) 

452(.90) 

452(.90) 

445(.89) 

458(.92) 

446(.89) 

too  low 

10(.02) 

28(.06) 

21(.04) 

36(.07) 

23(.05) 

43(.09) 

7.0 

too  high 

37(.07) 

25(.05) 

27(.05) 

18(.04) 

21(.04) 

10(.02) 

cover 

438(.88) 

434(.87) 

454(.91) 

442(.88) 

455(.91) 

439(.88) 

too  low 

25(.05) 

41(.08) 

19(.04) 

40(.08) 

24(.05) 

51(.10) 

9.0 

too  high 

47(.09) 

27(.05) 

26(.05) 

16(.03) 

23(.05) 

8(.02) 

cover 

357(.71) 

377(.75) 

455(.91) 

438(.88) 

454(.91) 

435(.87) 

too  low 

96(.19) 

96(.19) 

19(.04) 

46(.09) 

23(.05) 

57(.ll) 

11.0 

too  high 

59(.12) 

29(.06) 

26(.05) 

12(.02) 

23(.05) 

6(.01) 

cover 

233(.47) 

263(.53) 

455(.91) 

434(.87) 

453(.91) 

424(.85) 

too  low 

208(.42) 

208(.42) 

19(.04) 

54(.ll) 

24(.05) 

70(.14) 

13.0 

too  high 

59(.12) 

38(.OS) 

26(.05) 

11(02) 

23(.05) 

5(.01) 

cover 

154(.31) 

175(.35) 

454(.91) 

433(.87) 

453(.91) 

419(.84) 

too  low 

287(.57) 

287(.57) 

20(.04) 

56(11) 

24(.05) 

76(.15) 

15.0 

too  high 

83(.17) 

35(.07) 

26(.05) 

7(.01) 

23(.05) 

3(.01) 

cover 

73(.15) 

121(.24) 

456(.91) 

436(.87) 

453(.91) 

418(.84) 

too  low 

344(.89) 

344(.89) 

18(.04) 

57(.ll) 

24(.05) 

79(.16) 
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Table  31.     AVERAGE  AND  STANDARD  DEVIATION  OF  HALF  LENGTH(90% 
C.I,  MOD  1,  C=  1/10) 


TIME 

K.M.E 

M.L.E 

A.R.E 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

1.0 

.1121   .1105 

.0718   .0714 

.1276   .1187 

(.0064)  (.0062) 

(.0061)  (.0060) 

(.0812)  (.0704) 

3.0 

.1223   .1185 

.0972   .0953 

.1431   .1222 

(.0062)  (.0066) 

(.0042)  (.0044) 

(.0762)  (.0457) 

5.0 

.1123   .1046 

.0871   .0837 

.1300   .1021 

(.0132)  (.0142) 

(.0101)  (.0103) 

(.0765)  (.0379) 

7.0 

.0988   .0847 

.0688   .0644 

.1077   .0761 

(.0285)  (.0267) 

(.0134)  (.0135) 

(.0760)  (.0325) 

9.0 

.0806   .0640 

.0512   .0464 

.0861   .0533 

(.0437)  (.0363) 

(.0139)  (.0137) 

(.0749)  (.0268) 

11.0 

.0573   .0432 

.0369   .0321 

.0681   .0360 

(.0513)  (.0398) 

(.0128)  (.0122) 

(.0733)  (.0211) 

13.0 

.0416   .0305 

.0261   .0217 

.0540   .0237 

(.0506)  (.0380) 

(.0109)  (.0101) 

(.0716)  (.0160) 

15.0 

.0306   .0223 

.0182   .0145 

.0433   .0154 

(.0472)  (.0350) 

(.0089)  (.0079) 

(.0701)  (.0118) 
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Table  32.     TWO-SIDED  80  %  COVERAGE  FRACTION(MOD  1,  C=  1/10) 


TIME 

COV- 
ERAGE 

K.M.E 

M.L.E 

A.R.E 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

1.0 

too  high 

61(.12) 

47(.09) 

56(.ll) 

48(.10) 

35(.07) 

28(.06) 

cover 

404(.81) 

402(.80) 

400(.80) 

403(.81) 

396(.79) 

392(.78) 

too  low 

35(.07) 

51(.10) 

44(.09) 

49(.10) 

69(.14) 

80(.16) 

3.0 

too  high 

49(.10) 

38(.08) 

52(.10) 

42(.08) 

49(.10) 

37(.07) 

cover 

415(.83) 

417(.83) 

401(.80) 

403(.81) 

408(.82) 

399(.80) 

too  low 

36(.07) 

45(.09) 

47(.09) 

55(.ll) 

43(.09) 

64(.13) 

5.0 

too  high 

56(.ll) 

47(.09) 

52(.10) 

40(.08) 

49(.10) 

39(.08) 

cover 

411(.82) 

392(.78) 

397(.79) 

395(.79) 

409(.82) 

390(.78) 

too  low 

33(.07) 

61(.12) 

51(.10) 

65(.13) 

42(.08) 

71(.14) 

7.0 

too  high 

55(1 1) 

43(.09) 

51(.10) 

39(.08) 

47(.09) 

34(.07) 

cover 

410(.82) 

390(.78) 

395(.79) 

396(.79) 

407(.81) 

392(.78) 

too  low 

35(.07) 

67(.13) 

54(.ll) 

65(.13) 

46(.09) 

74(.15) 

9.0 

too  high 

71(.14) 

53(.ll) 

51(.10) 

34(.07) 

4S(.10) 

31(.06) 

cover 

333(.67) 

351(.70) 

394(.79) 

399(.80) 

407(.81) 

384(.77) 

too  low 

96(.19) 

96(.19) 

55(11) 

67(.13) 

45(.09) 

85(.17) 

11.0 

too  high 

79(.16) 

59(.12) 

53(.ll) 

31(.06) 

47(.09) 

25(.05) 

cover 

213(.43) 

233(.47) 

392(.78) 

400(.80) 

407(.81) 

381(.76) 

too  low 

208(.42) 

208(.42) 

55(.ll) 

69(.14) 

46(.09) 

94(.19) 

13.0 

too  high 

84(.17) 

56(.ll) 

55(.ll) 

29(.06) 

47(.09) 

24(.05) 

cover 

129(.26) 

157(.31) 

391(.78) 

398(.80) 

406(.81) 

379(.76) 

too  low 

287(.57) 

287(.57) 

54(.ll) 

73(.15) 

47(.09) 

97(.19) 

15.0 

too  high 

110(.22) 

58(.12) 

54(.ll) 

26(.05) 

48(.10) 

21  (.04) 

cover 

46(.09) 

98(.20) 

392(.78) 

398(.80) 

405(.81) 

373(.75) 

too  low 

344(.69) 

344(.69) 

54(.ll) 

76(.15) 

47(.09) 

106(.21) 
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Table  33.     AVERAGE  AND  STANDARD  DEVIATION  OF  HALF  LENGTH(80% 
C.I,  MOD  1,  C=  1/10) 


TIME 

K.M.E 

M.L.E 

A.R.E 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

1.0 

.0872   .0864 

.0559   .0557 

.0981   .0908 

(.0049)  (.0048) 

(.0047)  (.0047) 

(.0665)  (.0563) 

3.0 

.0945   .0927 

.0754   .0745 

.1070   .0929 

(.0050)  (.0051) 

(.0033)  (.0035) 

(.0593)  (.0357) 

5.0 

.0855   .0819 

.0670   .0654 

.0952  .0776 

(.0106)  (.0112) 

(.0080)  (.0081) 

(.0574)  (.0293) 

7.0 

.0729   .0664 

.0524  .0503 

.0771   .0578 

(.0216)  (.0209) 

(.0105)  (.0106) 

(.0571)  (.0250) 

9.0 

.0578   .0502 

.0385   .0362 

.0600   .0404 

(.0317)  (.0285) 

(.0108)  (.0107) 

(.0561)  (.0205) 

11.0 

.0403   .0339 

.0273   .0251 

.0461   .0273 

(.0364)  (.0312) 

(.0098)  (.0095) 

(.0548)  (.0161) 

13.0 

.0290   .0240 

.0190   .0169 

.0353   .0180 

(.0355)  (.0298) 

(.0083)  (.0079) 

(.0535)  (.0122) 

15.0 

.0213   .0175 

.0130   .0113 

.0274   .0116 

(.0330)  (.0275) 

(.0067)  (.0062) 

(.0524)  (.0089) 
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